Finite temperature pion vector form factors in 
Chiral Perturbation Theory 
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_-. - , Abstract 

We discuss the thermal behaviour of the pion vector form factors and calculate 
them in one-loop Chiral Perturbation Theory. The perturbative result is used to 
analyze the T-dependent electromagnetic pion charge radius, obtaining a rough 
estimate of the deconfinement critical temperature. Imposing thermal unitarity, we 
Q , generate the p resonance pole for the form factor in the center of mass frame. The p 

| peak height in the modulus of the form factor decreases for increasing temperature, 

while its width increases and the peak position is slightly shifted downwards for 
T ~ 150 MeV. These results point in the direction suggested by many analysis of 
^r) \ dilepton production data in relativistic heavy ion collisions. 

! PACS: 11.10. Wx, 12.39.Fe, 11.30.Rd, 25.75.-q., 12.38.Mh. 
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1 Introduction 

On going experiments on relativistic heavy ion collisions have attracted much attention 
over the past years. After the expected Quark-Gluon-Plasma cools down and hadronic 
matter forms, a correct description of the system involves QCD at temperatures below 

x ■ 

the chiral phase transition. In this regime, the medium constituents are predominantly 
light mesons, whose very low-energy dynamics is described by Chiral Perturbation The- 
ory (ChPT) [1, 2], the low-energy effective theory of QCD based on spontaneous chiral 
symmetry breaking. ChPT is a low-energy expansion performed in terms of p 2 /A^ with p 
a typical meson energy or temperature and A x ~ 1 GeV. It has been successfully applied 
to light meson dynamics and also to describe the low-T pion gas [3]. Recently, we have 
calculated in this framework the irn thermal scattering amplitude [4] and its unitarization 
[5] showing that chiral symmetry and unitarity alone provide a reasonable description of 
the thermal p and a resonances, without including them as explicit degrees of freedom 
(see also [6]). 

One salient observable in heavy ion collisions is the dilepton spectrum: e + e~ pairs 
are direct probes of the plasma evolution, since they do not interact from the production 
point to the detector. The observed spectrum in recent experiments performed at the SPS 
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collider for invariant masses between 200 MeV and 1 GeV differs significantly from vacuum 
hadronic emission models, showing a global enhancement for center of mass energies 
between 0.2-1 GeV [7]. In addition, the DLS collaboration at the BEVALAC accelerator 
has reported results with lighter ions and lower colliding energies per nucleon where the 
low-mass enhancement is even bigger and defies theoretical explanations [8]. 



The explanation of the dilepton SPS enhancement has been the subject of intense 
theoretical work [9, 10, 11, 12, 13]. This effect is particularly visible near the p mass, 
around which the spectrum flattens and is compatible with a widening of the p [9, 11, 
12, 14]. Ignoring baryon density effects, the main contribution to the dilepton spectrum 
at low energy stems from the annihilation of two thermal pions via the emission of a 
virtual photon 7r + 7r _ — > 7* — > e + e _ [15, 16] where the p is produced as an intermediate 
resonant state. Thus, the production rate for dileptons is governed by the pion form 
factor [15, 16] whose in-medium modification in the pion gas is key to the spectrum 
[9, 10, 11, 17]. Prior calculations have relied on model dependent input, the closest to 
our approach being the finite temperature form factor analysis in [17], a chiral model 
with resonances explicitly included as independent fields in the Vector Meson Dominance 
(VMD) framework. Another approach [18] uses a sum rule calculation to extrapolate to 
low energies from the perturbative results of QCD. In this work we first present a model 
independent study of the finite temperature effects on the pion form factor using ChPT. 
In particular, after a brief analysis of the different thermal form factors in Section 2, we 
obtain the one-loop ChPT thermal calculation, we study the temperature evolution of 
the pion electromagnetic charge and radius, obtaining a merely qualitative estimate of 
the critical temperature of deconfinement, and we check thermal unitarity. Second, in 
Section 4, we implement unitarity to describe the effect of the thermal variation of the p 
resonance mass and width in the form factor. In section 5 we summarize our main results. 
Comparing with our previous works [4, 5], we remark that, although the phase of the form 
factor is related with the amplitudes via unitarity, the modulus, which we calculate here 
in ChPT, is an independent object needed to describe relevant physical quantities such 
as the charge radius or the dilepton rate. 
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2 Finite temperature vector form factors 



At T 7^ 0, all physical quantities may depend on the fluid four-velocity, so that, with 
the usual choice of the fluid rest frame, Lorentz covariance is lost while spatial rotation 
covariance is still preserved. Therefore, the most general expression for the timelike vector 
form factor of charged pion or electromagnetic 7r + 7r~7 vertex, is 1 : 

(n + (p)n-(p')\V (0)\0} = q F t (S ,\S\,q ) 

(7r + (p)7T-(p')l^(0)|0> = q k F s (S ,\S\,q ) + S k q G s (S , \S\,q ) (1) 

with V M the electromagnetic current, S = p + p', q = p — p' and F t , F s , G s even functions 
in go (charge conjugation invariance). Note that the effect of Lorentz covariance breaking 
is twofold: on the one hand, the time and spatial components of the current may depend 
on the three different functions F t , F s and G s and, on the other hand, those functions may 
depend on three independent rotationally invariant variables, instead of just one invariant 
variable (s = S 2 ) as in the T = case. Terms containing ekijS l q j , allowed by rotation 
invariance and C are forbidden by parity. 

Of course, gauge invariance further restricts the above expressions, imposing relations 
between the functions F tj8 , G. When taking the divergence of the current, it must be also 
taken into account that the in-medium pion dispersion relation may be different from the 
vacuum one, i.e, the pole of the pion propagator is at p 2 = m 2 + g(po, \p\',T) with g a 
complex function. Altogether, this relates the form factors and the function g. However, 
we will content ourselves here with one-loop ChPT, which gives the real (and positive) 
leading 0(p 4 ) to g, that depends on T but not on energy and momentum [3]. That is, to 
0(p A ) the two pions in (1) can be treated as free with a T-dependent mass shift. Therefore, 
to 0(p A ) we have q • S — 0, that combined with gauge invariance (7T7t|<9 m ^|0) = in (1) 
yields: 

So - = |£| 2 G?> (2) 

where the superscript (1) means the NLO (one- loop) form factors. Remember that the 
(T-independent) leading order is just F^ = F^ = 1, G^ = 0. 

As we will see below, in one-loop ChPT G^ ^ and therefore there are two inde- 
pendent thermal form factors F^ ^ F^. Moreover, there is an additional simplification 

1 An equivalent parametrization is used in [17], where F s = F, G s = G' /{po — qo) and F t = F + G'(po + 
9o)/(po - 9o) + Gj (po - g ). 
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valid to one loop: the F t ,F s ,G s functions do not depend on q . The gauge invariance 
condition (2) coincides with the analysis performed in [17]. 



In this section we follow similar steps as in our pion scattering calculation [4]. We cal- 
culate the time-ordered product of two pion and one current fields in the imaginary time 
formalism (ITF) of Thermal Field Theory [19]. This we continue analytically for contin- 
uous external pion energies and connect to the form factor through the LSZ reduction 
formula. Such thermal amplitudes correspond to retarded real-time Green functions and 
have the correct analytic and unitarity structure [4, 5], properties of special interest in 
this study. 

In the ITF, the lagrangian and electromagnetic current coincide with those in ordinary 
T = ChPT. The thermal modifications arise upon replacing all zeroth momentum com- 
ponents by discrete frequencies k° — > iu n = 2ninT and the loop integrals by Matsubara 
sums, i.e., f^-—* iTJ2"^°- We draw in Figure 1 the diagrams (also the same as for 
T = 0) contributing to the form factors at NLO within the 577(2) chiral lagrangian [2]. 



Figure 1: Diagrams contributing to the pion form factors to one loop in ChPT. 

Diagram (a) gives, among other contributions, the imaginary part needed for unitarity 
(see section 3.2) while diagram (b) is proportional to the T-dependent tadpole which is 
real and affects the electromagnetic vertex through wavefunction renormalization. Their 
ITF contributions to the current expectation value, before performing the analytic con- 
tinuation are given by: 



3 One loop ChPT calculation 




(a) 





(3) 
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where So, k G 2iriTZ and D is the space-time dimension. As for diagram (c), the black 
box stands for the 0(p 4 ) lagrangian tree level contribution, and therefore T-independent, 
proportional to the low-energy constant l 6 [2]. 

The integrals appearing in (3) can be related to the T-dependent one-loop integrals 
discussed in the Appendix of [4]. We adhere to the notation and conventions of that 
work. Using the formulae there quoted and after analytic continuation, the form factors 
can be written in terms of three independent thermal integrals: J 2 , Jo, that are energy 
and momentum dependent, and the constant tadpole integral F, 



2. 



^\So,\S\) 



AFP(S ,\S\) = 
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(4) 



As in [4], for a T-dependent quantity we denote AH(T) = H(T) - H(0). Recall that 
all the D — > 4 UV divergences are contained in the T = part [2] that is finite and scale 
independent once expressed in terms of the finite and scale independent constant 1% and 
the subtracted T = loop integral J(s) = Jq(s) — Jq(0) 



F t (s) = F s (s) = l + J(sy^^ + ik - V3)S 



(5) 



6/2 96tt 2 /2 

We remark that our finite temperature additions (4) are written in Minkowski space- 
time for continuous So G R. In addition, we have used the on-shell condition p 2 = {p') 2 , 
valid to one loop. As we mentioned in section 2, the one-loop form factors do not depend 
on g - For instance, there are terms in (Vb) proportional to q ■ S = qoSo to this order 
that would otherwise give a nontrivial qo dependence. Finally, note that in the above 
expression and to this order we can use either the physical f n , m n , their tree level values 
or their T-dependent ones, since the differences are of higher order. We have chosen 
to write down our expressions in terms of the T = physical values /„- ~ 92.4 MeV, 
m T ~ 139.6 MeV. 



2 We have used Ji(S ,\S\) 



^Jo(S , \S\), that holds for any S , \S\. 



It is easy to check the consistency of our explicit one-loop form factors (4) with the 
gauge identity (2), showing how G^p ^ for arbitrary So, \ S\. Another interesting check 
concerning the unitarity of (4) in the center of mass frame will be analyzed in section 3.2. 

3.1 Pion electromagnetic charge and radius at T ^ 0. 

A direct prediction of our ChPT calculation is the pion electromagnetic static charge 
density at T ^ and low energies. Recall that at T = 0, the total pion charge and 
charge radius are (Q) = F(0) and (r 2 ) = 6F'(0)/F(0) where F(s) is the form factor. 
Thus, from (5), (Q) = 1 and (r 2 ) = (l e — 1)/I6ir 2 f%. In fact, this is the simplest way to 
estimate the value of T 6 in ChPT to one loop [2] giving l 6 = 16.5 ±1.1. Recent evaluations, 
including two loop corrections [20], obtain l e = 16.0 ± 0.5 ± 0.7, where the last error is 
purely theoretical and dominates the uncertainty. Concerning the charge radius, the latest 
experimental average quoted in the PDG is (r 2 ) = 0.45 ± 0.01 [21]. 

Likewise, at T ^ taking into account the Lorentz covariance breaking: 



2x 6 dF t (0,\S\) 



(r 2 )T = 



Qt d\S\< 



(6) 

\s\=o 



where Qt = lim F t (0, \S\;T). Note that the charge density is defined through V and 

|SHo+ 

hence F t in (1) must be used. We consider the spacelike radius as would be measured in 
t-channel e — it scattering and therefore So must be set to zero before \S\ (static limit). 
Analogous timelike and magnetic moment radii could also be defined. From equation (4) 
we find: 

2E 2 - ml 



rn 



2 



Q^) T = Wo + ^ r^-^[3(, 2 -. 2 )(, 2 + 2. 2 )^T) 



127T 2 / 2 J m7v E±^W^ 



+ E(2E i + 2ml ~ E K) 



dE 



(7) 



where ub{x;T) = [exp(a;/T) — l]^ 1 is the Bose-Einstein distribution function. 

Note that the positive charge decreases with temperature due to 0(T 2 ) corrections in 
the s — > limit. This is similar to the electric charge Debye screening in QED [19, 22]. 
For the negatively charged ir~ , the whole eq.(7) changes sign. However, in absolute value 
both charges decrease by the same amount and the gas remains neutral. We have plotted 
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(r 2 )r/ (r 2 )o in Figure 2. The T-correction is almost negligible below 100 MeV, where it 
decreases very slightly. For higher T, the radius increases considerably, the dominant 
contribution coming from the Qt screening discussed above. Since they rely on ChPT 
alone, these results are model independent. Our analysis confirms earlier results for the 
charge radius based on extrapolated QCD sum rules [18]. 

We get a rough estimate of the deconfinement temperature T c [22] when the pion 

1K(T), 

/ rz 



electromagnetic volume equals the inverse pion density, i.e, (47r/3)(r ) T 

r cPk i 

where n T (T) =3 j7^n B (J \k\ 2 + ml; T). Taking just (r 2 ) gives T c ~ 265 MeV, 



(2tt)* 

which is clearly too high, as commented in [22]. Thus, the thermal increase of {r 2 ) T 
reduces T c to a more realistic value. With our above result we get T c ~ 200 MeV, which 
in fact is closer to the chiral restoration transition temperature as estimated within ChPT 
[23]. The uncertainties in / 6 amount to ± 4 MeV. This critical temperature is below the 
temperature where Qt = and the radius diverges, which would correspond to a vanishing 
pion thermal mass and is clearly beyond the validity range of the chiral expansion. As a 
matter of fact, these estimates have to be regarded as merely qualitative. 




200 



Figure 2: The electromagnetic pion charge radius at T ^ 0. 

It has also been shown [2, 24] that the value of l 6 is saturated by the p resonance contri- 
bution, that, together with VMD, implies l e ~ 9Gir 2 f 2 /M 2 [24], up to chiral logarithms. 
Therefore, the electromagnetic radius behaves as M~ 2 . Thus, at T ^ 0, our previous 
results, still model-independent, would suggest an almost constant (though slightly in- 
creasing) M p (T) for very low T and a clearly decreasing M P (T) for T> 100 MeV. This 
is indeed the behaviour that we found in [5] and is also confirmed by our analysis in 
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section 4. Such behaviour for M P (T) at very low T is a prediction of chiral symmetry and 
resonance saturation, at least in the chiral limit [25]. 



3.2 Thermal perturbative unitarity in the center of mass frame 

Let us reduce our previous expressions to the center of mass frame (c.o.m.), i.e, p = —p' 
in (1). This amounts to look only to back to back e + e~ pairs in the dilepton spectrum. 
The dilepton rate is particularly simple in that case [16] although the three-momentum 
distribution may be important when performing more realistic analysis [11]. 

Therefore, we take the limit \S\ — > + in our previous expressions. Noting that 
J 2 = (Sg/4) J + Fp/2 + 0{\S\ 2 ) [4] we see in (4) that F t (1) , F s (1) and \S\ 2 G { P have a finite 
c.o.m. limit, which is reassuring. In addition, in the c.o.m. frame {p — p'){p + p') = 
implies S k q G s = S k (2Sty /S )G, = 2p k (\S\ 2 G s /S ) = q k (\S\ 2 G s /S ) so that \S\ 2 G S /S 
can be reabsorbed in F s in (1). Finally, we find: 

lim AF t (1 \S ,\S\) = lim AF^(S , \S\) = \(S 2 - Ami) AJ (S , 0) + AAFp 
|S|-o+ |5Ho+ 6J4 L 

Jim \S\ 2 GW(S ,\S\) = (8) 

|S|-0+ 

Therefore, in the c.o.m. frame there is only one T-dependent form factor, which we 
will just call F(So] T) for simplicity. The tadpole contribution Fp is real, while Jo contains 
a nonzero imaginary part required by unitarity. In fact, following the same steps as in [4]: 

Im F W (E + *e; T) = a T (E) E [~ ^ l (9) 

with 

a T (E) = a(E 2 )[l + 2n B (E/2;T)] (10) 



for positive energies above the two-pion threshold E > 2m- K . Here, a(s) = \Jl — 4m 2 /s 
is the two-pion phase space factor. 

At T = 0, (9) is the perturbative version of the form factor unitarity relation Im F(s) = 
a(s)F*(s)a u (s) where a 11 is the I = J = 1 partial wave projection of the irn scattering 
amplitude (to lowest order F — 1 and a n (s) = (s— 4m 2 )/9Qirf 2 ). At T 0, the correction 
factor 1 + 2n B = (I + n B ){l + ris) — n B riB is interpreted as enhanced phase space [26] 
due to the difference between induced emission and absorption processes [4]. Therefore, 
we find that the form factor satisfies a perturbative unitarity relation analogous to that 
in the T = case, but now in terms of a thermal phase space factor. The same happened 
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with the thermal amplitude in [4]. This is not only a consistency check of our calculation, 
but it will be the basis of our unitarization method used in the next section in order to 
generate the p pole in the form factor. 



4 Unitarization and applications 

The ChPT perturbative form factors analyzed in the previous sections provide the pre- 
diction of chiral symmetry to next to leading order at finite temperature in a model- 
independent way. However, by construction, they cannot reproduce a pole or resonant 
behaviour, which in the I = J = 1 channel corresponds to the p. Our approach here 
will be to construct a nonperturbative thermal form factor F imposing exact thermal 
unitarity with respect to the nonperturbative pion scattering amplitude. The latter is 
well approximated by the Inverse Amplitude Method (IAM) [27] derived at T ^ in [5] 
a IAM {E;T) = a 2 {E 2 )/[a 2 {E 2 ) - a 4 (E; T)] where a 2 and a 4 are the ChPT partial waves 
to tree and one-loop level respectively (we are suppressing the 11 superscript) calculated 
in [4]. The IAM amplitude is also exactly unitary, i.e, lma IAM = <7T\a IAM \ 2 so that our 
only physical input, apart from chiral symmetry, will be unitarity. It should be borne 
in mind that our approach of demanding exact thermal elastic unitarity is meant to be 
valid for energies and temperatures such that ub{E/2\T) remains small [5]. For the 
dilepton spectrum this means that for typical freeze-out temperatures T ~ 150 MeV, 
our approach is rather accurate around the p scale, which is precisely where unitarity is 
saturated. Near the pion pair threshold, corrections to pion propagation not included in 
this work (nominally of 0(p e ) in ChPT) have been conjectured to be more important 
[10, 11, 16, 17]. 

If we take F oc a IAM with a real proportionality constant we readily guarantee the 
exact unitarity condition ImF = axF*a IAM and also that both the amplitude and the 
form factor have the same poles in the complex plane as well as the same complex phase 
(phase shift), as it happens for T = 0. Imposing the correct low energy perturbative 
ChPT expansion F — 1 + F^ + ■ ■ ■ fixes the real proportionality constant, so that: 

F(F .T)- 1 + **FW(E;T) a 2 (E 2 ) 
1 ' ) a 2 (E 2 )+Rea 4 (E;T) a 2 (E 2 )-a 4 (E;T)- 1 ) 

valid for E > 2m n , where perturbative unitarity holds. The above formula was developed 
at T = in the strongly interacting electroweak spontaneous symmetry breaking sector 
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[28]. At low energies it reproduces the chiral expansion up to plus terms of higher 
order, which should be smaller. Thus we expect that it should reproduce also the low 
energy data with an l 6 value slightly different, but reasonably close, to that used in 
the previous section with the pure one-loop ChPT. All that remains is to adjust the 
undetermined low-energy constant Iq to zero-temperature experimental data and the finite 
T behaviour follows as a prediction. For a IAM and hence for a<i and 04 we are using the very 
same calculation given in [4, 5] but note that in order to obtain the unitarized \F(E;T)\ 
from (11) we need the real part of the form factor, whose perturbative calculation we have 
carried out in section 3. In Figure 3a we show how the p data [29] are nicely described 
by the resulting phase shift in the 11 channel, common to the amplitude and the form 
factor. This fixes the form factor phase. In Figure 3b. we plot \F(E;0)\ 2 for 0.3 GeV 
< E < 1 GeV compared with the data in [30]. The solid line has been obtained with 
Iq = 18, that describes data reasonably well and is compatible with the perturbative value 
quoted in section 3.1. For illustration we also provide curves with Iq = 17 (dotted) and 
Iq = 19 (dashed). Observe that the resonant behaviour around M p ~ 770 MeV is clearly 
reproduced in both figures. In fact, note that the unitarized formula gives exactly the 
form factor phase and therefore the position and width of the p pole is very accurately 
reproduced. However, the modulus is more subject to perturbative uncertainties, as 
reflected in Figure 3. 




300 400 500 600 700 800 900 1000 300 400 500 600 700 800 900 1000 



E (MeV) E (MeV) 

Figure 3: a) We show the phase of the a IAM for elastic amplitude of irn scattering in the p 
channel. The data comes from [29] b) The unitarized form factor at temperatures T = 0, 
for Iq = 18 (solid line), Iq = 17 (dotted line) and Iq = 19 (dashed line). The data comes 
from [30]. 

Next, in Figure 4 we have plotted \F(E;T)\ 2 for different temperatures. We observe 
that the form factor decreases and widens with temperature. The mass position of the 



10 



peak moves slightly to the right and then drops in the T = 150 MeV curve, consistently 
with our discussion on the pion electromagnetic radius in section 3.1. Note that, by 
construction, the form factor (11) has a peak exactly at the same place and with the 
same width as the amplitude, that we had already studied in [5], and where the p pole 
moved further away from the axis, explaining the strong flattening of the form factor in 
this work. 

Our results broadly agree with [17]. Note however that we have not introduced explicit 
resonances, the physical assumptions are just chiral symmetry and unitarity. Our T = 
peak also falls a little bit short of the data but is much closer than [17]. Since we only 
deal with two pions, and therefore only the p resonance, this could be partially due to uj 
contamination in the data coming from e + e~ annihilation (first reference in [30]). Indeed, 
the lowest data point at the peak, closer to our curves, comes from r decay (second 
reference in [30]), where no uj can be produced. The two-loop calculation [31, 20] can also 
improve the situation, but this is beyond our scope at T ^ 0. 



300 400 500 600 700 800 900 1000 
E (MeV) 

Figure 4: The unitarized form factor at temperatures T = 0, 100, 150 MeV. 

The modulus of the form factor enters directly in the dilepton rate from pion anni- 
hilation. The simplest approach is to use kinetic theory to leading order, so that the 
dilepton rate is written as an integral containing the pion distribution functions and 
\F\ 2 [15, 16]. In the c.o.m. and in thermal equilibrium the rate is simply proportional to 
n 2 3 (E/2; T)\F(E; T)\ 2 [16]. However, the equilibrium result is not enough to fit the exper- 
imental data. The dynamics of the plasma expansion, the contribution of other channels 
and the experimental acceptance have to be properly accounted for. Let us remark that, 
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as pointed out in [10, 17], it might be crucial that there is a quark-hadron phase transition 
and a rather long-lived mixed phase in order that medium modification effects coming 
from pion annihilation are sizable in the dilepton yield. In [11] it has been found that 
to fit the dilepton data reasonably well, medium effects need to be included in the form 
factor. Otherwise, the theoretical prediction would exceed the dilepton data around M p 
which means that medium effects should decrease the height of the p peak. Moreover, the 
form factor is expected to spread by a factor of two at T ~ 150 MeV [11] and the peak 
position to shift to lower energies, as in earlier treatments [9] . More elaborated space-time 
analysis including also baryon effects [12] conclude that the p spreading allows to explain 
the CERES data, without a sizable mass decrease. All these effects are qualitatively visi- 
ble in our theoretical result for the form factor, our peak mass position shift being smaller 
than assumed for instance in [11]. In this respect, it is worth pointing out that a recent 
measurement by the STAR collaboration [32] of the p° in-medium properties directly in 
p° — > 7r + 7r~ shows a softer p-mass decrease (of about -40 MeV) and therefore much closer 
to the claims in [12] and also to the size of our predictions here. This is an important 
measurement, as chiral symmetry restoration requires that p and a\ become degenerate 
[33] so that this predicts a sharp a x mass decrease. Further lowering of the p mass as 
measured in pion observables can be achieved theoretically due to various minor medium 
effects [34] such as collisional broadening, the Boltzmann factor and rescattering at later 
stages of the collision (lower temperatures). 

5 Conclusions 

We have analyzed the pion vector form factors at finite temperature in Chiral Perturbation 
Theory. The general structure of the form factors at T ^ allows in principle for three 
different form factors. However, the gauge Ward-Takahashi identity constrains them, 
relating the form factors to the in-medium pion dispersion relation. 

Our explicit one-loop calculation gives the two different form factors not tied by gauge 
invariance with the correct T = limit. In the center of mass frame, the two form factors 
coincide and satisfy a perturbative thermal unitarity relation in terms of a thermal phase 
space, consistently with our previous results on the thermal tttt scattering amplitude. 
Using only ChPT, we have also studied how the effective charge and charge radius of 
the pion change with temperature. The effective charge is screened with T, while the 
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radius is almost constant for low T and then increases. Up to here, these results rely only 
on thermal field theory and ChPT and are therefore model independent. We have also 
checked that this behaviour is consistent with the expected p mass thermal behaviour. 
Our form factor has also allowed us to obtain a naive estimate of the deconfinement 
temperature. 

On a second stage, by imposing exact thermal unitarity while respecting the ChPT low 
energy expansion, we construct a nonperturbative thermal form factor that reproduces 
previous theoretical analysis and whose behaviour is qualitatively compatible with the 
observed dilepton spectrum. At the typical freeze-out temperatures of T ~ 150 MeV, 
our result predicts a peak height decrease and a spread of the form factor around the p 
region, as expected from dilepton data. In addition, the position of the peak is slightly 
shifted to lower mass. We have arrived to our result imposing only chiral symmetry and 
thermal unitarity as physical assumptions. 

The unitarization method discussed here is limited to the center of mass frame. For 
future work, it would be interesting to extend these ideas including the three-momentum 
dependence needed for the dilepton analysis. According to our discussion above, this 
would need to account simultaneously for the effect of different form factors and the 
dispersion relation. Still, an analysis of back-to-back lepton data by ongoing experimental 
collaborations would allow comparison with these very simple and powerful theoretical 
results. In addition, a more realistic study should also include the space-time evolution 
of the plasma and baryon density effects. Work is in progress along these directions. 
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